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Abstract. We present a wide grid of models for the structure and transmission properties of warm absorbers in 
active galactic nuclei (AGN). Contrary to commonly used constant density models, our absorbing cloud is assumed 
to be under constant total (gas plus radiation) pressure. This assumption implies the coexistence of material at 
different temperatures and ionization states, which is a natural consequence of pressure and thermal equilibrium. 
Our photoionization code allows us to compute the profiles of the density, the temperature, the gas pressure, the 
radiation pressure and the ionization state across the cloud, and to calculate the radiative transfer of continuum 
and lines including Compton scattering. Therefore, equivalent widths of both saturated and unsaturated lines 
are properly modeled. For each pair of the incident spectrum slope and the ionization parameter at the cloud 
surface there is a natural upper limit to the total column densities of the cloud due to thermal instabilities. These 
maximum values are comparable to the observational constraints on the column density of warm absorbers which 
may give support to constant total pressure models. In all models we note considerable absorption around 6.4 keV 
which modifies the intrinsic relativistically broadened iron line profile originating in an accretion disk illuminated 
atmosphere. Our models can be applied to fitting the spectroscopic data from the XMM-Newton and Chandra 
satellites. 
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1. Introduction 

The first X-ray absorption feature due to ionized heavy 
elements was recognized by Halpern (1984) in the X-ray 
spectrum of the Syl galaxy MR 2251-178 observed by 
the EINSTEIN satellite. The author related the observed 
jump in flux around 1 keV to the absorption edge of Oviii. 
Therefore, as Halpern concluded, X-rays emitted from the 
central region of an active galaxy traveling toward an ob¬ 
server encounter a “warm absorber” — material with an 
electron temperature lower than the temperature of colli- 
sionally ionized gas with a similar level of ionization. 

Results from EXOSAT, ROSAT, GINGA, ASGA, and 
Beppo-SAX satellites showed that warm absorbers were 
common among Seyfert galaxies. Nandra & Pounds (1994) 
have suggested that more than 50% of Syls contain a 
warm absorber, and later results confirmed such occur¬ 
rence rate (e.g. Reynolds, 1997; George at al. 1998). Those 
conclusions were based exclusively on the detection of ab¬ 
sorption edges (for instance in MCG-6-30-15, Nandra & 
Pounds 1992). The detector areas and the spectral resolu¬ 
tion of those satellites did not allow us to see any absorp¬ 
tion lines from highly ionized species. 

Send offprint requests to: A. Rozahska (agata@camk.edu.pl) 


The situation changed since 1999, when large X-ray 
telescopes Ghandra and XMM-Newton started to operate, 
with their X-ray grating instruments working in the en¬ 
ergy range up to almost 10 keV. Several tens of absorption 
or emission lines were observed and identified in NGG 3783 
(Kaspi el al. 2002, Behar el al. 2003, Netzer et al. 2003, 
Krongold el al. 2003), NGG 5548 (Kaastra el al. 2002), 
NGC 1068 (Kinkhabwala et al. 2002), NGC 7469 (Blustin 
et al. 2003), MGG-6-30-15 (Turner et al. 2004). In other 
objects the number of fitted lines is lower and/or their de¬ 
tection is less firm, but the results still strongly support 
the presence of a warm absorber in those Seyferts (NGG 
4051, Collinge et al. 2001; Mrk 509, Yaqoob el al. 2003; 
TonS180, Rozahska el al. 2004). Warm absorber lines were 
even detected in a distant blazar at z = 4.4 (Worsley et 
al. 2004). 

Detected lines are basically consistent with the uni¬ 
fication scheme of AGN based on the presence of the 
dusty/molecular torus. According to that scheme, Syl are 
objects seen face-on and Sy2 are objects seen edge-on 
(Antonucci & Miller 1985; for application to X-ray band 
see Mushotzky, Done & Pounds 1993). The major charac¬ 
teristic of the X-ray spectra of Sy2 galaxies is that they are 
strongly absorbed at low energies, and that their emission 
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lines, especially the iron Ka line, have large equivalent 
widths (Turner et al. 1997) since their intensity is mea¬ 
sured with respect to the heavily obscured direct contin¬ 
uum (Weaver & Reynolds 1998) or only with respect to the 
scattered continuum (in Compton-thick objects), as dis¬ 
cussed by Bassani et al. (1999). In some cases of Compton- 
thin Sy2 obscuration may be due to the host galaxy and 
unrelated to disk/torus orientation (Matt 2000, Guainazzi 
et al. 2001) The spectra of Syl galaxies are predomi¬ 
nantly featureless, show little low energy absorption, and 
the equivalent width of the iron Ka line is much lower 
than in Sy2. Some Syl show relativistically broadened 
iron Ka line profiles which additionally supports the view 
that in Syl we directly observe the innermost part of the 
nucleus. The scattering medium may be identical to the 
warm absorber and we have side view of this medium in 
Sy2 galaxies while we see the nucleus through it in Syl. 
Indeed, X-ray spectra of some Sy2 galaxies show many 
narrow emission lines, (e.g. NGC 1068, Kinkhabwala et 
al. 2002; NGC 4507, Matt et al. 2004; Mkn 3, Pounds & 
Page 2005; NGC 4151, Schurch et al. 2004) indicating that 
ionized material may extend beyond the shielding torus. 
In Syl spectra mostly absorption lines are observed (e.g. 
NGC 3783, Kaspi et al. 2002: NGC 5548, Kaastra et al. 
2002; Mrk 509, Yaqoob et al. 2003) without strong ab¬ 
sorption of continuum, suggesting that in those objects 
we see directly the nucleus through the ionized plasma. 
Some contribution of emission lines in Syl X-ray spec¬ 
tra is also expected (e.g. Netzer 1993, Collin, Dumont & 
Godet 2004) and seen in the data (e.g. NGC 3783, Behar 
et al. 2003; NGC 7469, Scott et al. 2005). However, the 
constraints from Sy2 galaxies for the ionized medium can¬ 
not be used directly to Syl galaxies since in Sy2 galaxies 
we may observe only the outer part of the plasma distri¬ 
bution, due to the obscuration by the torus, while warm 
absorber features may come predominantly from the inner 
part. 

A majority of these absorption lines shows a veloc¬ 
ity shift of the order of a few hundreds km/s (Kaspi 
et al. 2001, Kaastra et al. 2002) suggesting that the 
warm absorber is outflowing. A strong high velocity out¬ 
flow was recently reported in several objects, mostly 
radio-quiet quasars (e.g. 60000 km/s and 120000 km/s 
in APM 08279-f5255, Chartas et al. 2002; 23000 km/s 
in PG1211-I-143, Pounds et al. 2003a; 63000 km/s in 
PG0844-t349, Pounds et al. 2003b, 30000 km/s in IRAS 
I3I97-I627, Dadina & Cappi 2004; 26,000 km/s in PG 
I404-I-226, Dasgupta et al. 2005). However, the results are 
based on a few detected lines, so line identifications and, 
consequently, the determined high outflow velocities can 
be questioned (Kaspi 2004a). The origin of the outflow and 
its geometry is still under discussion (Grenshaw, Kraemer 
& George 2003a, Blustin et al. 2005). 

It is more difficult to estimate the radial distance of 
the warm absorber from the nucleus. Absorption lines 
most probably form somewhere between the broad line 
region (BLR) and the narrow line region (NLR) (i.e. from 
about 0.01-0.1 up to 10 pc from nucleus; see Grenshaw et 


al. 2003a, Blustin et al. 2005). Variability in the overall 
warm absorber properties was reported for a few sources 
(MR 2251-178, Halpern 1984, Kaspi et al. 2004b; MGG 
-6-30-15, Reynolds et al. 1995; H1419-I-480, Barcons et al. 
2003; NGC 4395, Shih et al. 2003; NGC 3516, Netzer et 
al. 2002). Netzer et al. (2002) and Barcons et al. (2003) 
concluded that the changes observed are consistent with 
varying ionization of the gas. The lack of short-timescale 
(days) response of the warm absorber to the change of the 
continuum can be used to put lower limits to the warm 
absorber distance (e.g. 0.5 - 2.8 pc for NGC 3783, Behar et 
al. 2003; similar limits were given by Netzer et al. 2003). 
Krongold et al. (2005) detected a response of the warm 
absorber to the changes of the continuum on a timescale 
of 31 days in NGC 3783, deriving an upper limit of 6 
pc for the distance of the warm absorber. The shortest 
variability timescale of ~ 10^ s has been detected in the 
warm absorber in MCG -6-30-15 (Otani et al. 1996; see 
also Turner et al. 2004) locating the plasma responsible 
for the Oviii edge within the distance of 10^^ cm from the 
nucleus. The spectral analysis of this source indicates that 
the highly ionized warm absorber is dust-free, with dust 
contribution in this source coming from a distant zone, 
hundreds of pc from the nucleus (Ballantyne et al. 2003). 

The same medium is most probably responsible for 
narrow absorption lines seen in the UV spectra of many 
AGN (for a review, see Crenshaw et al. 2003a). Some kine¬ 
matic components discovered in UV coincide with those 
discovered in soft X-rays but for other components no such 
correspondence is seen (e.g. Behar et al. 2003 for NGC 
3783, Crenshaw et al. 2003b for NGC 5548, Kaspi et al. 
2004b for MR 2251-178; Scott et al. 2005 for NGC 7469; 
Gabel et al. 2005 for NGG 3783). However, the resolving 
power of the UV observations is high {R ~ 20000) while 
X-ray data resolution is much lower {R ~ 1000), which 
makes the comparison difficult, as discussed by Crenshaw 
et al. (2003a). Analysis usually suggests that UV and X- 
ray absorption features are consistent with arising in the 
same gas, but with stratified ionization (e.g. Barcons et 
al. 2003, Kaspi et al. 2004b, Scott et al. 2005, Gabel et al. 
2005). Since the absorption features appear in the profiles 
of the broad emission lines like CiV, and are occasionally 
deep, this serves as an argument that the absorbing region 
is located outside the BLR. 

The column density of the warm absorber is generally 
estimated to be about cm“^, and absorbing gas 

contains heavy elements mostly in the form of helium- 
and hydrogen-like ions. However, accurate measurements 
of the column density are quite complex. Most determi¬ 
nations are based on detection of absorption edges, but in 
some cases edges are undetectable while lines are clearly 
seen (Kaastra el al. 2002, Rozahska et al. 2004). Also the 
estimates of column densities from edges do not always 
confirm estimates derived from the absorption line analy¬ 
sis (Kaspi et al. 2002). 

The ionization state of the gas required to reproduce 
the observed absorption or emission lines seems to be 
quite complex. In many objects predictions based on a 
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single cloud at one specific ionization state cannot explain 
the data, so an absorbing material is modeled using at 
least two photoionization regions, which are required to 
explain presence of lines from the matter in different ion¬ 
ization states (see in the case of Syl: NGC 4051, Collinge 
et al. 2001; NGC 5548, Kaastra et al. 2002; NGC 3783, 
Kaspi et al. 2002, Netzer et al. 2003 and Krongold et al. 
2005, H0557-385, Ashton et al. 2005). The same conclu¬ 
sion was drawn from fitting of the absorbed continuum 
(strongly absorbed Syl: Mkn 304, Piconcelli et al. 2004; IC 
4329A, Steenbrugge et al. 2005a,b; in the case of a dwarf 
galaxy with an active nucleus: NGC 4395, Shih, Iwasawa 
& Fabian 2003; and in case of Sy2: NGC 4507, Matt et al. 
2004). 

The idea of the warm absorber being under the 
constant pressure was developed even before Chandra 
and XMM-Newton satellites were lunched (Netzer 1993; 
Krolik & Kriss 1995; Netzer 1996). It is well known that 
cold/warm material irradiated by hard X-rays should be 
strongly stratified, and eventual thermal instabilities lead 
to its clumping (McKee Krolik & Tarter 1981). If ther¬ 
mal instabilities are strong, none of the existing photoion¬ 
ization codes can describe unstable zone. This is because 
radiative transfer codes are based on unique density and 
temperature profiles, and do not accept situations where 
for one value of optical depth the solution gives three dif¬ 
ferent values of temperatures and densities. 

Therefore many models of the warm absorber being 
under constant pressure do assume that those two dis¬ 
crete phases already exist. Usually two or three zones at 
different constant densities are assumed to have the same 
dynamical ionization parameter, S, which is the ratio of 
ionization pressure to the gas pressure (McKee Krolik & 
Tarter 1981). In such situation, calculations of transfer 
of X-ray radiation through two constant density zones are 
done separately and then spectrum is merged together de¬ 
pending on covering factor (Netzer 1993). 

In this paper we solve the situation where X-ray ra¬ 
diation is no so hard and strong, that separation on two 
phases takes place. Instead, our warm absorber is strongly 
stratified and radiation passes through different densities 
and ionization stages. Theoretically this situation was con¬ 
sidered by Krolik & Kriss (2001), Krolik (2002). The strat¬ 
ification of our warm absorber is determined physically by 
radiative properties, and the only assumption which we 
make is constant total pressure within a cloud. We use 
the photoionization code titan developed by Dumont et 
al. (2000) (see Dumont et al. 2003 for implementation of 
Accelerated Lambda Iteration method) to compute the 
synthetic spectra for a systematic set of model parame¬ 
ters. The advantage of our calculations is that we compute 
the full radiative transfer of continuum and lines taking 
care also of lines which are saturated. Recent observations 
suggest that saturated lines are often present in warm ab¬ 
sorber (Kaspi et al. 2002, Rozahska et al.2004). The trans¬ 
fer is done in the stratified medium, and the density profile 
is determined self-consistently with the radiation transfer 
to fulfill the condition of the pressure equilibrium across 


the cloud. We aim understanding these AGN which show 
clear absorption lines in their Chandra or XMM spectra 
so the warm absorber is located in the line of sight to the 
observer. We concentrate on modeling the transmission 
spectra and absorption lines therefore our models can be 
used for Syl galaxies. 

We describe our model in Sec. |2 and El while results 
are presented in Sec. □ □ and El Sec. 13 describes com¬ 
parison of our models to observations and Sec. El contains 
conclusion remarks. 

2. Description of the model 

The main assumption of our model is that the structure 
of a single warm absorber cloud is determined by the con¬ 
dition of constant total, Pgas + Prad, pressure. The value 
of this pressure is self-consistently determined by solving 
non-LTE ionization equilibrium, thermal equilibrium and 
radiative transfer throughout the cloud thus providing a 
density and ionization profile inside a cloud. 

We model an irradiated cloud using a plane parallel ap¬ 
proximation. The total column density of the cloud, Nh, 
is a free parameter in our model. 

The second model parameter is the ionization param¬ 
eter, determined at the cloud surface 

L AttF 

= -> ( 1 ) 

noA no 

where L is the luminosity of the central source, no is the 
number density at the cloud surface, at its illuminated 
side, and R is the distance of the cloud from the central 
source. There is no difference if we use small ^ or dy¬ 
namical big S ionization parameter which is defined as 
S = F/{cPgas) since both parameters are used only to 
specify the amount of X-ray flux illuminating the cloud 
surface. In our models both ionization parameters change 
with the optical depth due to the change in the density 
and in the irradiating flux, and this stratification is com¬ 
puted self consistently. 

The X-ray illuminating continuum is assumed to have 
a power-law shape, as it is usually observed in many 
Seyferts, and the power-law photon index, F, is a model 
parameter. 

The photoionization code titan developed by Dumont 
et al. (2000) solves the full radiative transfer of the ion¬ 
izing continuum in two stream approximation. The trans¬ 
fer method for the lines and continuum used in this code 
has been modified to use the accelerated lambda iteration 
(ALI) method, which is especially well suited to rapid con¬ 
vergence in the cores of lines which may be optically thick 
(Dumont et al. 2003). The computations are done assum¬ 
ing complete redistribution function in the lines. Partial 
redistribution is mimicked by a Doppler profile for some 
of the most intense resonant lines. 

The comparison of proper ALI computations with es¬ 
cape probability approximation was shown in the recent 
papers by Dumont et al. (2003) and Collin, Dumont, 
Godet (2004). In the case of emission lines the errors in 
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the line fluxes are typically of the order of 30% for col¬ 
umn densities 10^° cm“^ and a factor of five for column 
densities of 10^^ cm“^, compared to calculations using the 
escape probability formalism. 

TITAN includes all relevant physical processes from 
each ion level. The population of each level is computed 
solving the set of ionization equations coupled with the 
set of statistical equations describing the excitation equi¬ 
librium. However, some low ionization states are more 
roughly treated. 

The following processes are taken into account: radia¬ 
tive and collisional ionization in all levels, recombinations 
in all levels, and radiative and collisional excitations and 
de-excitations for all transitions. 

We consider ten main heavy elements including iron. 
TITAN transfers “only” 900 lines while some photoioniza¬ 
tion codes (XSTAR, cloudy) collect over 4000 lines or 
more. On the other hand, titan is also designed to cal¬ 
culate radiative transfer in Compton thick media, which 
basically means that we are able to compute transfer 
of optically thick lines through any stratified cloud, and 
also through an atmosphere in hydrostatic equilibrium 
(Rozahska et al. 2002). The same code is applicable for 
a whole range of column densities (from small to large), 
and a whole range of ionization states. 

The Compton heating/cooling balance through the 
cloud is taken into account since for hard photon index, 
r = 1.5, Compton heating is comparable to photoioniza¬ 
tion heating. Compton balance of photons up to 26 keV 
is computed in titan while Compton balance of photons 
between 26 up to 100 keV is computed using code NOAR 
(Dumont, Abrassart, & Collin 2000). The NOAR is a ra¬ 
diative transfer code, based on the Monte Carlo method, 
to compute a Compton heating profile for a plane strati¬ 
fied atmosphere. This Compton heating/cooling profile is 
read by the titan and taken into consideration for solv¬ 
ing the radiative transfer. Usually three iterations between 
NOAR and TITAN are enough to determine Compton heat¬ 
ing/cooling curve. For the sets of our models with the 
same F the Compton heating does not vary much with ^ 
and of the warm absorber. Hence, we only derived 

a few representative Compton profiles, for different F and 
we adopted them to all other models to perform our com¬ 
putations. Comptonization of line photons is included. 

Generally our code models the reflected and outward 
spectrum from the slab with assumption that both are 
emitted isotropically. In Fig.^we present all spectra avail¬ 
able for one particular cloud. The advantage of our code 
is that we do not separate absorbed and emitted spectra 
which are derived from two different parts of the solu¬ 
tion of radiative transfer. This procedure is done in many 
previous attempts starting from Netzer (1993). In titan 
we solve transfer of absorption and emission continuum 
and lines simultaneously producing outward spectrum (see 
Fig.Hl), which is a mixture of absorption and emission fea¬ 
tures. Nevertheless, if a warm absorber cloud is on the line 
of sight to the observer, we see outward cloud spectrum 



Fig. 1. All spectra available from one considered cloud 
(^ = 10®, Ntot = 10^^). Incident power-law is indicated 
by dotted line and is assumed to be isotropic. Reflected 
spectrum is emitted isotropically on the same side of the 
cloud, while outward spectrum is emitted isotropically on 
the opposite side of the cloud. Transmitted spectrum is 
outward spectrum in normal direction to the observer. 

emitted only in the direction normal to the slab, i.e. point¬ 
ing to the observer. 

In this first paper we want to examine only the absorp¬ 
tion features, therefore the final spectra presented in the 
section below are purely absorption spectra, i.e., outward 
in the direction toward the observer (transmitted). This is 
equivalent to the situation where the covering factor of the 
warm absorber is small, but a single warm absorber cloud 
is large and located on the line of sight to the observer. 
The model is applicable to Seyfert 1 galaxies since in those 
objects we observe mostly absorption lines. Although the 
absorbed spectra do not contain any emission lines, the 
temperature and ionization structure were computed with 
full transfer of absorbed, emitted and reflected radiation. 
Reflected and outward spectra in other directions will be 
treated more carefully in the next paper. 

In this paper generally we are discussing a grid of mod¬ 
els with zero velocities to see the main trends in the prop¬ 
erties of the absorbing matter. But we also present two 
models with non-zero turbulent velocity the same for all 
ions to trace how it affects final equivalent widths of lines 
(see also a special case considered in Rozahska et al. 2004). 

3. Grid of parameters and numerical limitation 

All computations are done assuming that the surface of 
a warm absorber has a number density uq = 10^^ cm“®. 
The exact value of this parameter seems to be unimpor¬ 
tant. We have checked that adopting a lower value of the 
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surface density does not change the cloud structure. The 
ionization parameter ^ gives a direct information on the 
flux of the ionizing radiation which hits the surface of the 
warm absorber. In our model we vary ^ from 10^ up to 
10® erg cm s“^. 

We assume that the ionizing continuum is coming from 
the center of an active nucleus and has a power law shape. 
The photon index varies from 1.5 up to 2.5 as observed for 
the hard X-ray band. The illuminating X-ray continuum 
in all models extends from 0.01 keV up to 100 keV. 

When computing warm absorber models under con¬ 
stant pressure we deal with the well known problem of an 
illuminated atmosphere with heavy elements. If the inci¬ 
dent spectrum is hard enough and/or the abundance of 
heavy elements is solar or higher, a thermal instability de¬ 
velops (Fields 1965, Krolik, Mckee & Tarter 1981) which 
affects the conditions in the surface layers of the irradiated 
stellar or accretion disk atmospheres, including AGN disks 
(e.g. Rozahska & Czerny 1996). 

The instability occurs when for a single value of the 
dynamical ionization parameter S = Ftot/cPgas the en¬ 
ergy balance equation has three solutions with different 
temperatures and densities. Therefore, in a constant den¬ 
sity medium it is not possible to obtain any instability but 
they develop if the density of the medium can adjust to 
the local conditions. Among these three solutions, one is 
unstable, and two others are stable, but there is no stable 
solution which is continuous through the whole irradiated 
layer. In our case we assume constant pressure instead of 
constant density and we may expect thermal instabilities 
in the clouds with total column density high enough for 
the temperature to drop considerably inside the cloud (see 
Sec. EJ. 

Generally, for each set of F and ^ we consider a range 
of values for the logarithm of the column density, log 
starting with 21. Whenever possible, we extended our grid 
to log Nh — 23.5. However, for specific values of the initial 
ionization parameter and/or steepness of X-ray spectrum 
we were not able to achieve thermal equilibrium. There 
is a maximum total column density, for which in¬ 

stabilities are so strong that computations failed. For this 
reason our grid of total column densities is not the same 
for each set of F and These maximum values of the total 
column density are shown in Fig.[21for all combinations of 
F and 

The maximum value rises with an increase of 

ntot 

the ionization parameter but the relation is not strictly 
monotonic, and the dependence on the spectral slope is 
even more complex. The most difficult case was F = 1.5 
and ^ = 10^, where a temperature fall appeared even for 
low total column densities, and we were able to find a 
stable cloud only up to the N}itot = 3.16 x 10^^ cm“^. 
Under those special conditions, line cooling due to a single 
iron ion strongly dominated all other cooling mechanisms, 
and the medium was optically thick in the emitted lines. 
Therefore, the local temperature strongly depended on the 
presence or absence of a layer behind the considered zone. 
The cooling of the outer zone served as a heating mecha- 



log f 

Fig. 2. The maximum total column density versus ioniza¬ 
tion parameter for each photon index. For a given value 
of F and ^ the structure of clouds with column densities 
higher than cannot be determined because of ther¬ 

mal instabilities. 

nism to the inner zone, compensating exactly the losses. 
As a result, the Thomson depth of the hot zone depended 
roughly linearly on the assumed ■ 

In this paper we also study the case when absorbing 
matter has turbulent velocity. For this purpose we con¬ 
sider two cases, one with Vturb — 100 km/s, and one with 
Vturb = 300 km/s, the same for each ion. 

4. Structure of the warm absorber 

We expect that at the illuminated face of a cloud the tem¬ 
perature will be high and density low, and a high tempera¬ 
ture equilibrium on the “S curve” (Krolik, Mckee & Tarter 
1981) will be reached. This layer is expected to be optically 
thin despite some lines which may be saturated. Going 
deeper inside the warm absorber the illuminating contin¬ 
uum will be absorbed by heavy elements and we expect 
that the temperature structure will show a strong decrease 
with accompanying increase of the density, while the ion¬ 
ization state decreases, but not so rapidly. Therefore after 
this geometrically extended optically thin layer, we will 
get geometrically thin layers with strong temperature and 
density gradients becoming optically thicker. 

The temperature structures of all computed clouds are 
given in Fig. 0 Models with F = 1.5 are presented in 
panels a, b, c for ^ = 10®, 10"^ and 10® respectively. The 
same sequence of but for F = 2, is shown in panels d, 
e, f, and the softest intrinsic spectrum with photon index 
2.5 is presented in g, h and i panels. The density structure 
for exactly the same set of models is presented in Fig. 0] 

For each set of F and ^ we present results for a grid of 
column densities starting from to about 

^Htot ’ depending on the appearance of thermal instabili¬ 
ties. We use a 0.5 step size in logarithm of the total column 
density. 

In each panel (except for b, where F = 1.5 and ^ = 10"^) 
we see the same trend in the structure while increasing 
the value of the total column density. For the lowest to¬ 
tal column densities the temperature does not decrease 
strongly, because the cloud is thin enough to exist only 
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Fig. 3. The temperature profiles within the irradiated 
clouds at constant pressure. The different panels show 
the dependence on the total column density, starting with 
logNHt^t = 21, with step sizes of 0.5. The step size de¬ 
creases to 0.1 when approaching the value of . 

Panels a, b, c show the results for ^ = 10®, 10^ and 10®, 
and for T = 1.5. Panels d, e ,f show respectively the same 
ionization parameters but for P = 2. The last three panels, 
g, h, i, present the cases for P = 2.5. 
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Fig. 4. The density profiles within the irradiated clouds 
at constant pressure. The different panels show the de¬ 
pendence on the total column density, starting with 
log Nntot =21, with step sizes of 0.5. The step size de¬ 
creases to 0.1 when approaching the value of logA^^““. 
Panels a, b, c show the results for ^ = 10®, 10^ and 10®, 
and for P = 1.5. Panels d, e, f show respectively the same 
ionization parameters but for P = 2. The last three panels, 
g, h, i, present the cases for P = 2.5. 


in the high temperature equilibrium. In such a case the 
warm absorber is well modeled by a constant density slab 
with low density of the order of n = 10®® cm“®. A temper¬ 
ature drop appears when Nutot increases, and the highly 
ionized hot layer is complemented by geometrically thin, 
dense zones reaching even n = 5 x 10®'® cm“®. Since the 
transition to solutions with the outer cool cloud layer oc¬ 
curs rapidly, we had to decrease the step in logNHtat 
0.1, for the last two curves in almost each panel in order 
to resolve correctly the structure. 

The exceptional case ofP = 1.5, ^ = 10® (panel b 
in Fig. 01 which was discussed at the end of Section 01 
shows a strong temperature jump even for Njitot = 
cm“^. The dense layer on the back of the illuminated face 
has an unphysical temperature bump. Such a layer has a 
non negligible optical thickness and it is seen in modeled 
spectra, which possesses strong absorption (see Section 
below). 

In Fig. 01 we present thermal equilibrium curves, T(S), 
for several clouds. The value of dynamical ionization pa¬ 
rameter across the cloud was computed using relation 
S = Prad/Pgas ■ From this figure it is clearly seen the differ¬ 
ence of our models with those when two constant density 
slabs are used to fit a data. We propose warm absorber as 
a single cloud which passes through all ionization states, 


and each layer interacts radiatively with another one as it 
is usually done in classical atmospheric calculations. 

5. Modeled spectra 

The transmitted spectra for the whole set of models are 
presented in Figs. E117I and IHl for F = 1.5,2 and 2.5 re¬ 
spectively. All figures are given in the range of energies 
between 0.1 and 4 keV, where most lines are present. The 
energy ranges of iron lines from 6 keV up to 8 keV are 
presented and discussed in Sec. El 

For high ionization parameters and low total column 
densities the spectra are featureless since the matter is al¬ 
most fully ionized and only changes in the continuum are 
caused by Comptonization. This effect is clearly seen in 
the upper panels of Figs. El and 0 where the whole outgo¬ 
ing continuum is lowered. The Compton heating reaches 
almost 50% of total heating, therefore the determination 
of the exact value of Compton heating was very important 
in such cases and we did it iteratively using the Monte 
Carlo code NOAR. 

In the spectra corresponding to higher total column 
densities we observe many absorption features. The spec¬ 
trum in the energy range from 0.8 keV up to 2 keV is 
especially strongly absorbed due to many transitions of 
highly ionized heavy elements. 
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Fig. 5. The thermal equilibrium curves for the clouds at 
constant pressure. Clouds with T = 1.5 are presented in 
the left panel, with T = 2.0 in the middle panel, and with 
r = 2.5 in the right panel. Lines have following meaning: 
dotted line - ^ = 10®, dashed line - ^ = 10"^, and solid line 
for clouds with ^ = 10^. In all cases we show clouds with 
maximum total column density. 



1 

V (kcV) 


Fig. 6. Set of transmission spectra for the photon index 
of the incident continuum L = 1.5 and ^ = 10® (upper 
panel), ^ = 10“^ (middle panel) and ^ = 10® (lower panel). 
The results for a range of XogNnt^t are shown in each 
panel starting from log Nh = 21. up to logNf^^J^. 


The comparison of spectra for the same total column 
density of the cloud but for different ionization parameters 
is shown in Fig.|3 The upper panel presents the results for 
slabs with NHtot — 10^^ cm“^ illuminated by radiation of 
photon index T = 1.5. The different spectra from the top 
to the bottom represent the cases for ^ = 10®, 10"'^ and 10® 



1 

V (keV) 


Fig. 7. Set of transmission spectra for the photon index 
of the incident continuum T = 2.0 and ^ = 10® (upper 
panel), ^ = 10^ (middle panel) and ^ = 10® (lower panel). 
The results for a range of logA^/z^^j are shown in each 
panel starting from logNn = 21. up to logA^^““. 


respectively. The same sequence of ionization parameters 
is presented in the middle panel but for Nutot = 10^^ 
cm“^ and T = 2.0, and in the lower panel but for = 

3.16 X 10^^ cm“^ and T = 2.5. 

The spectrum in the case of T = 1.5 and ^ = 10^ is 
again exceptional because of the numerical problems men¬ 
tioned above. In this spectrum many absorption edges are 
seen. This is because the optically thick layer on the back 
of illuminated cloud is also geometrically thick and the 
absorption is much stronger. Additionally, the highly ion¬ 
ized optically thin layer on the front of illuminated cloud 
is thin enough to not interact with X-ray radiation at all, 
explaining the lack of lines above 1 keV in this spectrum. 

Other spectra show an increase in the number of lines 
with decreasing ionization parameter. Also less energetic 
lines (below 1 keV) are much more numerous for lower 
This is because the back side of an illuminated cloud 
reaches lower temperatures and higher densities, and also 
lower ionization states for lower 

5.1. Absorption edges 

The detection of absorption edges in the spectra of AGN 
puts direct constraints to the total column density of the 
warm absorber. Therefore, we show in Fig.^|the strength 
of the main ionization edges in our synthetic spectra for 
comparison with observations. Models with T = 1.5 are 
presented in panels a, b, c for ^ = 10®, 10“^ and 10® respec¬ 
tively. The same sequence of but for T = 2 is shown 
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Fig. 8. Set of transmission spectra for the photon index 
of the incident continuum F = 2.5 and ^ = 10® (upper 
panel), ^ = 10“^ (middle panel) and ^ = 10® (lower panel). 
The results for a range of logNnt^t are shown in each 
panel starting from log Nh = 21. up to log 


in panels d, e, f, and the softest intrinsic spectrum with 
photon index 2.5 is presented in g, h and i panels. 

The strength of an edge is measured by its total op¬ 
tical thickness defined by total absorption plus scatter¬ 
ing opacity coefficient multiplied by the density and inte¬ 
grated over the geometrical size of cloud. We show only 
Ovii (filled triangles), Oviii (open squares), Cv (stars), 
and Cvi (open circles) edges above Tedge = 0 . 003 , which 
can be resolved using presently working X-ray satellites. 

Again panel b shows the exceptional case, and we see 
that the ionization edges are strong for the reason men¬ 
tioned in the previous Section. In the other cases the edges 
become visible only when the total column densities are 
higher than 10 ^® cm“^ for f = 10 ® (panels a, d, g). This 
minimum total column density decreases with decreasing 
ionization parameter reaching the value of 10^^ cm“^ for 
^ = 10 ®. 

5.2. Strength of the absorption lines 

Several physical conditions influence the amount of energy 
absorbed in lines: the heavy element abundance, tempera¬ 
ture of the medium, the hardness of the illuminating con¬ 
tinuum, the column densities of particular ions, and the 
turbulent velocity of the absorbing matter. It is not ob¬ 
vious which condition is the most important. The curve 
of growth analysis should be done, but at present in our 
basic models we didn’t allow for any velocity dispersion in 
the absorbing gas. 



1 

V (keV) 


Fig. 9. Comparison of spectra for the same Nutot^ but for 
different ionization parameters. In upper panel we com¬ 
pare clouds of Nutat = 10^^ cm“^ illuminated by radiation 
with photon index T = 1.5, middle panel Nutot = 10^® 
cm“^ and F = 2.0, and lower panel = 3.16 x 10^^ 

and F = 2.5. In each panel lines from the top to the bot¬ 
tom present cases for ^ = 10 ®, 10 ^ and 10 ® respectively. 


In order to study the strength of absorption lines in our 
spectra, in Figures cniniini and d we show expanded 
views of the spectra corresponding to F = 2 (see Fig. 13. 
The range between 0.3 and 0.7 keV is shown in Fig. 1111 
between 0.7 and 1.5 keV in Fig. d from 1.5 up to 3.5 keV 
is presented in Fig. d and finally the range from 6.2 up 
to 9 keV is shown in Figd 

All lines are weak (see TableCJ, with equivalent widths 
smaller than 1 eV, and only in a few cases (Sixiv Tj/a, 
Sxvi Lya, Fexxvi Lya up to Ly^, also some fluorescent 
iron lines) do they reach a few eVs. For a decreasing ion¬ 
ization parameter f (going from upper to bottom panels 
on each figure) we see that the lines become stronger and 
are present even for lower total column densities. 

Some lines do not change very much with increasing 
column densities. These lines are predominantly from the 
most highly ionized species. Lines from ions of slightly 
lower ionization degree are not present in cases when the 
total column density is low, because the dense zone does 
not exist on the back of illuminated cloud. 

Several modeled lines do agree with those detected by 
X-ray satellites as is shown in Tab. ^ But most of them 
are slightly lower which may indicate that the observed 
fines are blended. We do expect this effect with present 
resolution of X-ray gratings. Also, it may indicate that 
material is in turbulent motion. 

Therefore, for one specific case (F = 2.0,^ = lO"*) we 
calculate the models with the maximum value of Af^“( 





















A. Rozariska R. Goosmann A.-M. Dumont B. Czerny: Modeling the Warm Absorber in AGN 9 

Table 1. Equivalent widths in [eV] of absorption lines from clouds with different ionization parameters at the surface. 


r 

2.0 

2.5 

2.0 

2.0 

Observed 

ref. 

C 

10® 

10'‘ 

10® 

10® 

10-* 

10^ 



log(iVfftoti) 

23.5 

23.5 

23.5 

23.2 

23.5 

23. 








Vturb =100 km/s 

Vturb =300 km/s 



Cvi Lya 

0.18 

0.24 

0.25 

0.31 

0.5 

1.27 

0.98 ±0.33 

(2) 

Nvil Lya 

0.19 

0.28 

0.31 

0.4 

0.63 

1.55 

1.33 ±0.24 

(2) 

Ovill Lya 

0.54 

0.65 

0.79 

0.98 

1.13 

2.97 

1.85 ±0.55 

(1) 

Oviii (0.7743) 

0.38 

0.44 

0.44 

0.38 

1.05 

2.76 

1.77 ±0.44 

(1) 

Fexvii (0.8211) 

0.002 

0.35 

0.93 

0.90 

0.7 

1.24 

1.42 ±0.28 

(1) 

Fexix (0.8266) 

0.002 

0.33 

0.52 

0.47 

0.97 

2.50 

0.83 ±0.24 

(1) 

Fexviii (0.8731) 

0.01 

0.45 

0.88 

0.47 

1.0 

2.37 

1.52 ±0.34 

(1) 

Fexvii (0.8984) 

- 

0.1 

0.23 

0.29 

0.2 

0.2 

1.43 ±0.52 

(2) 

Fexix (0.9184) 

0.1 

0.96 

1.17 

0.96 

1.47 

3.48 

1.84 ± 0.54 blend 

(2) 

Neix (0.9218) 

0.03 

0.33 

0.63 

1.25 

0.8 

1.48 

4.13 ±0.48 

(1) 








2.6 ± 1.48 

(2) 








i.oo_Q 3Q 

(3) 








i.uy 0.34 

(4) 

Fexx (0.953) 

0.2 

1.22 

1.22 

0.97 

1.69 

3.95 

4.88 ± 0.42 blend 

(1) 

Nex Lya 

0.6 

0.77 

0.86 

0.78 

1.54 

4.08 

3.20 ± 0.24 blend 

(1) 

Neix (1.127) 

0.003 

0.13 

0.27 

0.27 

0.2 

0.2 

all complex 


Fexvii (1.127) 

- 

0.04 

0.2 

0.22 

0.07 

0.07 

4.66 ±0.33 

(1) 

Fexxiii (1.127) 

0.964 

3.82 

2.86 

1.12 

4.54 

9.55 

2.15 ±0.82 

(2) 

Nex (1.209) 

0.36 

0.54 

0.55 

0.46 

1.40 

3.42 

2.3 ±0.24 

(1) 

Mgxil Lya 

0.74 

0.99 

1.08 

0.82 

2.09 

5.45 

4.39 ±0.24 

(1) 

Fexxiv (1.556) 

0.43 

0.5 

0.44 

0.1 

1.69 

3.94 

1.84 ±0.23 

(1) 

Mgxii (1.745) 

0.4 

0.67 

0.7 

0.56 

1.83 

4.09 

1.75 ±0.22 

(1) 

Sixiii (1.865) 

0.46 

1.47 

1.92 

1.91 

2.66 

6.6 

4.18 ±0.20 

(1) 








3.37 ± 1.68 

(2) 








1 4Q+0.76 

-L-"±y o.73 

(4) 

Sixiv Lya 

1.26 

1.78 

1.80 

1.13 

3.12 

7.99 

6.61 ±0.26 

(1) 








3.87 ± 1.94 

(2) 









(4) 

Sixiv (2.370) 

0.8 

0.97 

0.94 

0.71 

2.72 

6.64 

5.03 ± 0.72 blend 

(1) 

Sxv (2.461) 

0.72 

2.23 

2.56 

2.02 

3.62 

8.73 

4.49 ± 0.59 

(1) 

Sixiv (2.500) 

0.48 

0.84 

0.85 

0.63 

2.33 

4.48 

1.96 ±0.71 

(1) 

Sxvi Lya (2.611) 

1.51 

1.93 

1.81 

0.96 

3.68 

9.3 

5.88 ±0.66 

(1) 








q Qr+2.86 

O.OO 2.80 

(4) 

Sxv (2.883) 

0.3 

0.93 

1.04 

0.86 

2.87 

6.15 

1.81 ±0.67 

(1) 

Sxv (3.032) 

0.13 


0.9 

0.72 

2.19 

3.31 

2.82 ±0.82 

(1) 

Sxvi (3.100) 

0.9 

1.11 

1.04 

0.69 

3.10 

6.48 

2.87 ± 0.78 blend 

(1) 

Fexxv (6.700) 

6.74 

9.0 

5.35 

0.25 

11.99 

25.19 

13.76 ±5.07 

(1) 

Fexxv (7.881) 

3.37 

3.57 

2.22 

0.047 

8.25 

17.75 

33.06 ± 19.04 

(1) 


(1) NGC 3783, 900 ksec Chandra Kaspi et al. 2002; (2) NGC 5548, 86.4 ksec Chandra Kaastra et al. 2002; (3) Mrk 509, 59 ksec 
Chandra Yaqoob et al. 2003; (4) NGC 4051, 81.5 ksec Chandra Collinge et al. 2001. 


assuming two values of the turbulent velocity: 100 km/s 
and 300 km/s. The turbulence affects the thermal struc¬ 
ture of the cloud (see Fig. [T^ . The interior of the turbu¬ 
lent cloud is hotter than the interior of a ’static’ cloud, 
and cloud interior is less dense. The continuum is slightly 
modified in the 2-3 keV range but absorption lines are 
naturally much more intense (see Fig. CEl). 


Quantitative results for the line intensities are shown 
in Tab. Q] We see that the line equivalent widths increase 
by a factor of a few. The same trend was already seen 
in a constant pressure models of Rozahska et al. (2004). 
Higher value of the turbulent velocity even generally leads 
to overprediction of the line strength. However, some pre- 
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log N„^ log log 


Fig. 10. The strength of absorption edges for irradiated 
clouds at constant pressure. The different panels show the 
dependence on the total column density, starting with 
logNHt^t = 21, with step size of 0.5. The step size de¬ 
creases to 0.1 when approaching the value 
Panels a, b, c show the results for ^ = 10®, 10"^ and 10®, 
and for T = 1.5. Panels d, e, f show respectively the same 
ionization parameters but for P = 2. Last three panels: g, 
h, i present the cases for P = 2.5. Filled triangles denote 
Ovii edges, open squares are Oviii, stars are Cv, and open 
circles are Cvi edges. 

dieted iron lines are still definitively weaker than in the 
data. 

We plan to consider more systematically the effects of 
turbulent velocities in our future work. Metal abundance 
is another free parameter in our model which should be 
studied more carefully in the next paper. 

6. Influence of the warm absorber on the 

relativistic iron line 

A broad, relativistically smeared iron line was discovered 
in a number of sources, MCG -6-30-15 being the first dis¬ 
covered and the best example of this phenomenon (Tanaka 
et al. 1995). It is produced by the reflection of the hard X- 
ray continuum on the disk surface very close to the black 
hole. If the radiation originating at the disk passes through 
the warm absorber before it reaches the observer the spec¬ 
tral shape of the line is modified. 

To show this effect, we considered an example of the 
broad disk line. We have assumed that the disk is illumi¬ 
nated by a random distribution of flares following Czerny 
et al. (2004). Each flare produces a hot spot below, and 
each spot emits locally like an illuminated disk atmo¬ 
sphere. The final emissivity of the disk was computed by 



0.3 0.4 0.5 0.6 

V (keV) 


Fig. 11. The expanded version of Fig. [7| (P = 2) in the 
softest energy range between 0.3 and 0.7 keV. Models for 
different ^ = 10®, 10'^, and 10® are presented in upper, 
middle and lower panel respectively. The lowest continu¬ 
ous lines in each panel correspond to appropriate 

for the set of P and ^ as seen if Fig. |21 



0.8 1 1.2 1.4 


V (keV) 

Fig. 12. The same as Fig.^Jbut for energy range between 
0.7 and 1.5 keV. 

integrating over all disk radii. The black hole was assumed 
to be rotating, with a Kerr parameter a = 0.95 and an in¬ 
clination to the observer of * = 30°. All effects of general 
relativity were included. 

The spectral shape of the reflected disk component 
from this model (without any effect of the warm absorber) 
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1.5 2 2.5 3 

v (keV) 


Fig. 13. The same as Fig. [TTlbut for energy range between 
1.5 and 3.5 keV. 



6.5 7 7.5 8 8.5 9 

1' (keV) 

Fig. 14. The same as Fig.^Jbut for energy range between 
6.2 and 9 keV. 


is shown in Fig.^l^s the dashed line in both a and b pan¬ 
els. Since the local reflection comes from the partially ion¬ 
ized disk surface, several components of the iron line com¬ 
plex are still visible in the resulting spectrum, although 
they are relativistically smeared. However, if we see the 
nucleus through the warm absorber with a total column 
density Nfjtot — 4-5 x 10^^ and an ionization parameter 
^ = 10"* (second spectrum from the bottom in FiEf. 1141 in 
the middle panel), the overall spectrum is suppressed and 



Fig. 15. The structure of warm absorber in case of tur¬ 
bulent velocity for F = 2.0, ^ = 10^ and maximum col¬ 
umn density Panel a) shows temperature struc¬ 

ture, b) thermal equilibrium curve T(S) and panel c) den¬ 
sity structure. Solid line represents the case with Vturb = 0, 
dotted line with Vturb = 100 km/s, and dashed line with 
Vturb = 300 km/s . 




6.5 7 7.5 8 8.5 9 

V (keV) 


Fig. 16. Set of transmission spectra in case of turbulent 
velocity for F = 2.0, ^ = 10^ and maximum column den¬ 
sity Solid line represents the case with Vturb = 0, 

dotted line with vturb = 100 km/s, and dashed line with 
Vturb = 300 km/s . 


the narrow absorption features appear superimposed onto 
the broad emission line profile. 

In panel a of Fig. El we show the result (solid line) as 
calculated with the resolution used in our computations 
E/A.E = 390 , while in b panel we show the same result 
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Fig. 17. An example of the effect of absorption onto the 
shape of a relativistically smeared disk line. The unab¬ 
sorbed disk spectrum is marked by dashed line, while the 
solid line represents spectrum with two resolutions: a res¬ 
olution of 390 in panel a, and the resolution of Chandra 
HEG at 6 keV i.e. 200 in panel b. 

degraded down to the resolution of E/AE = 200 charac¬ 
teristic of the Chandra HEG instrument at 6 keV. 

We see that a strong absorption modifies the overall 
shape of the broad iron line. The absorbed shape of line 
seems to be more symmetric and closer to a Gaussian 
shape. Also the fitted representative energy (mean energy 
of the line) would be lower than derived for the intrinsic 
spectrum. The situation is therefore similar to the case of 
UV spectra where we see narrow absorption features in 
BLR lines. 

7. Comparison with observations 

Recent observations of Syl or Syl.5 galaxies show the 
presence of absorption lines in soft X-ray spectra from 
Chandra or XMM satellites. The number of detected lines 
depends on the resolution and sensitivity of the detec¬ 
tor, and the best observations were obtained for NGC 
3783 (Kaspi el al. 2002, Behar el al. 2003, Netzer et al. 
2003, Krongold el al. 2003). A 900 ksec observation al¬ 
lowed for the detection of over few tens of lines. In the 
cases of about 90 ksec observations or shorter, only a 
few lines were detected and fitted like in case of Mrk 
509, TonSlSO, PG1211-bl43, or NGC 7469 (Blustin et al. 
2003). Therefore, it is hard to fit our models with those 
observations, but we make some comparison by examining 
parameters already fitted by other authors. 

First of all, the data suggest that the total column 
density of the warm absorber is never higher than 10^^ 
cm“^ (Piconcelli et al. 2004). In our models we have a 


physical limit of the column densities because of the pres¬ 
ence of thermal instabilities. For each spectral index and 
ionization parameter there is a maximum value of the col¬ 
umn density for which warm absorber exists (see Fig. El). 
In all cases the values of this maximum total column den¬ 
sities agree with the column densities determined from 
observations of different object. 

As we mentioned above, in a stratified medium it is 
difficult to determine the plasma parameters by compar¬ 
ing directly the EWs of individual lines. Therefore, in 
Tab. El we present the ratios of the EWs of lines for the 
same ion. In the first column names of ratios are given, 
in the second column those ratios are computed theoret¬ 
ically assuming that lines are optically thin (see formula 
2. in Rozahska et al. 2004). We have listed only those line 
ratios which were detected by X-ray satellites. Columns 
from 3 up to 7 present the observed ratios assuming the 
most pessimistic error, when the numerator has its maxi¬ 
mal value, and the denominator has its minimal value (see 

Tab.EJ. 

The observed line ratios are usually much lower than 
predicted in the optically thin medium limit. Therefore, 
observations strongly suggest that lines are optically thick 
and radiative transfer is required to predict equivalent 
widths. 

The same ratios computed for our models are pre¬ 
sented in Tab. El For each set of models with spectral 
index and ionization parameter, the ratios are of the or¬ 
der of those observed when we consider clouds with higher 
total column densities. We conclude here that when the 
warm absorber have a total column density higher than 
10^^ cm“^ the resonance lines of the main ions become 
optically thick and radiative transfer should be used to 
model them. For instance, the most popular Li/a lines of: 
Oviii, Nex, Mgxii, Sixiv, and Sxvi for ^ = 10"‘ F = 2.0 
and logIV_f/j^j=23.5 have optical depths integrated on the 
whole cloud equal respectively: 977, 363, 234, 436, and 
216. 

8. Discussion 

We have calculated a grid of models of warm absorber 
with density stratification determined by the condition of 
the constant pressure. The assumption of constant pres¬ 
sure instead of constant density was motivated by recent 
observations of AGN. In most observed objects the ion¬ 
ization states implied by the observed lines span a large 
range, which cannot be accounted for a single photoion- 
ized region, or by collisionally ionized matter. At least two 
photoionized absorbing regions are required to fit the data 
and the properties of these two regions are consistent with 
pressure equilibrium 

Our models represent an essential improvement in the 
description of a single cloud over the models presented in 
the literature. All previous models adopted constant den¬ 
sity approximation for the medium (Netzer 1993, Kaastra 
et al. 2002, Kinkhabwala et al. 2002, Krongold et al. 2005), 
or even constant temperature (Krolik & Kriss 1995). On 
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Table 2. Ratios of Equivalent Widths of lines from the same ion in the optically thin approximation. For comparison 
the results from observations are given if available. We calculate the error in the most pessimistic way, when the 
numerator has its maximal value EW+AEW, and the denominator has its minimal value EW-AEW. 


Name of Ratio 

Theory 

NGC 3783 

NGC 5548 

Mrk 509 

NGC 4051 

Ton S180 


fij/fik 

Kaspi et al. 

Kaastra et al. 

Yaqoob et al. 

Collinge et al. 

Rozahska et al. 



2002 

2002 

2003 

2001 

2004 

CV(o.3079)/(0.3545) 

4.6 

- 

3.85 ±23.15 

- 

- 

- 

CVl(Lya)/(Lyl3) 

8.49 

- 

2.33 ±2.99 

- 

- 

- 

CVl(Lya)/(Lyj) 

14.34 

- 

3.77 ± 14.94 

- 

- 

- 

CVl(Lya)/(LyS) 

29.86 

- 

- 

- 

- 


^yil(Lya)/(Lyf3) 

5.25 

1.53 ±4.39 

12.09 ± 1.18 

0.86 ± 0.98 

- 

- 

'^^^^(Lya)/(Lyy) 

14.34 

3.5 ±8.3 

- 

- 

- 

- 

^Vll(^Lya)/(LyS) 

29.92 

1.48 ±2.27 

- 

- 

- 

- 

OVII(o,574)/(o.6656) 

4.75 

0.95 ± 1.76 

1.32 ±0.58 

- 

1.11 ± 1.45 

- 

OVII(o,574)/(o.6977) 

1.26 

0.73 ± 1.01 

2.13 ±2.21 

- 

1.43 ± 2.48 

- 

OVlll(Lya)/ (Lyl3) 

5.26 

1.04 ±0.75 

1.16 ±0.55 

- 

2.06 ± 1.37 

- 

OVlll(^Lya) / (Ly^) 

14.34 

0.89 ±0.55 

1.54 ± 1.17 

- 

- 

- 

OVlll(^Lyo,)/ (LyS) 

29.92 

0.71 ±0.41 

7.36 ± 1.46 

- 

- 

- 

MgXI(i,352)/(l,579) 

4.84 

1.15 ±0.19 

- 

- 

- 

- 

MgXI(i,352)/(l,660) 

12.9 

1.46 ±0.34 

- 

- 

- 

- 

SiXIII(i,865)/(2,182) 

4.97 

1.73 ±0.51 

- 

- 

- 

- 

SiXIII(i,865)/(2.294) 

13.09 

2.19 ± 1.37 

- 

- 

- 

- 

SiXIV(7.j,Q,)/(7.j,^) 

14.34 

3.37 ±2.12 

- 

- 

- 

- 

SXV(2,461)/(2.883) 

4.76 

2.48 ± 1.98 

- 

- 

- 

- 

SXV(2,461)/(3.032) 

13.18 

1.58 ±0.95 

- 

- 

- 

- 

FeXVII(o,8211)/(0.8984) 

9.51 

0.47 ±0.24 

0.88 ±0.99 

- 

- 

- 

FeXVIpo,8211)/(1.008) 

3.01 

0.67 ±0.22 

- 

- 

- 

- 

FeXXV(6,700)/(7.881) 

4.95 

0.42 ± 0.93 

- 

- 

- 

- 


the other hand, we do not consider here a global picture 
including the dynamics of the warm absorber flow which 
form a separate and a major issue (see Chelouche & Netzer 
2005 and the references therein). 

We have calculated the full radiative transfer of the 
illuminating radiation in continuum and in lines through 
the plane-parallel density-stratified slabs of different to¬ 
tal column densities. For a low total column density the 
cloud is mostly ionized and optically thin, and does not 
differ from the constant density model. But for a higher 
total column density the optically thick dense zone arises 
at the back of the illuminated cloud, and the temperature 
falls dramatically. Therefore, absorbing matter contains 
zones of different ionization states coexisting under con¬ 
stant pressure. 

For each set of models with spectral index E and sur¬ 
face ionization parameter ^ there is a maximum total col¬ 
umn density for which the cloud is thermally stable. 

For higher column densities thermal instabilities do not al¬ 
low the calculations to converge. Interestingly, is of 

the order of the maximum column density derived from 
X-ray observations of different AGN. It gives additional 
support to the idea of the warm absorber clouds being in 
pressure equilibrium. 

The modeled lines are usually optically thick and their 
equivalent widths are of the order of the observed values. 
We conclude that the observed warm absorbers at column 
densities of the order of 10^^ cm“^ and higher possess 


saturated absorption lines and that full radiative transfer 
is required to model their equivalent widths properly. For 
full evaluation of this problem we will perform a curve of 
growth analysis in our future work. 

Generally lines are easier to detect than ionization 
edges and this tendency is observed in several AGN 
(Kaastra et al. 2002, Rozahska et al. 2004). Our models 
predict that edges are not observable up to column den¬ 
sities about 10^^ cm“^ and even higher, depending on E 
and 

The most interesting result of our computations is the 
shape of the spectrum around an iron Ka line. In almost 
all models there are strong and narrow absorption fea¬ 
tures due to highly ionized iron ions, above Fexvii. Such 
absorption affects the shape of the broad iron emission 
line possibly originating in illuminated disk atmosphere. 
For instance, relativistic broad iron line profile after pass¬ 
ing through warm absorber becomes disrupted into three 
narrower profiles, which can be fitted by Gaussian shape. 
Such lines were reported in several AGN as presented 
by Yaqoob & Padmanabhan (2004) (see also NGC 3783 
Reeves et al. 2004). 

In this paper we consider only single cloud but obser¬ 
vations suggest that absorbing material forms distribution 
of clouds with covering factor which may depend on ve¬ 
locity of clouds (de Kool, Korista & Arav 2002). However, 
theoretical modeling of cloud distribution with complex 
velocity field is complicated. Also, the study of the trans- 
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mission spectra, including absorption lines, do not give us 
any diagnostic of the matter along other directions than 
the line of sight. The titan code can be used, however, to 
study the radiative transfer in all directions. We address 
both issues to the future work. 
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Table 3. Ratios of Equivalent widths of lines from the same ion computed for our models. We have chosen two column densities for each set of T and ^ where 
more less transition from optically thin to optically thick cloud occurs . 


r 

1.5 

2.0 

2.5 

C 

10 

5 

10 

4 

10® 

10® 

10 

4 

10 

3 

10 

5 

10 

4 

10 

3 

logiNutot) 

23.7 

23.8 

21. 

22.5 

22. 

22.6 

21.5 

23.5 

21.5 

23.5 

21.5 

23.5 

21.5 

23.5 

21.5 

23.4 

21.5 

23.2 

CV(o,3079)/(0.3545) 

4.57 

1.12 

1.27 

3.88 

4.49 

3.32 

4.61 

4.61 

4.61 

4.56 

4.61 

1.34 

4.60 

1.72 

4.61 

2.15 

4.61 

2.40 

CVl(Lya)/(Lyl3) 

5.09 

2.04 

1.24 

1.55 

1.32 

1.40 

5.27 

2.54 

5.25 

1.15 

3.52 

1.16 

5.19 

1.19 

4.68 

1.16 

2.13 

2.03 

CYl(Lya)/(Ly-f) 

13.64 

2.23 

1.20 

2.03 

2.17 

1.67 

14.35 

6.06 

14.3 

1.35 

9.05 

1.31 

14.37 

1.39 

12.77 

1.56 

4.97 

2.45 

CVl(Lya)/(LyS) 

28.61 

2.30 

1.17 

2.62 

3.81 

2.11 

29.89 

12.12 

29.79 

1.69 

18.54 

1.52 

29.65 

1.61 

26.26 

1.57 

9.75 

2.80 

^Yll(Lya)/(Lyl3) 

5.25 

1.73 

0.97 

1.35 

1.64 

1.26 

5.26 

3.12 

5.27 

1.18 

4.08 

117 

5.23 

1.24 

4.87 

1.20 

2.58 

1.86 

(Lyy) 

13.91 

1.75 

1.09 

1.82 

3.23 

1.59 

14.36 

7.8 

14.34 

1.46 

10.76 

1.36 

14.64 

1.50 

13.50 

1.45 

6.40 

2.27 

^Vll(^Lya)/(LyS) 

29.61 

1.82 

1.61 

2.32 

6.25 

2.18 

29.89 

15.88 

28.57 

1.99 

22.18 

1.62 

29.79 

1.74 

27.38 

1.69 

12.57 

2.56 

OVII(o,574)/(o.6656) 

4.80 

3.55 

2.96 

2.95 

2.75 

3.02 

4.76 

4.63 

4.77 

1.76 

4.68 

2.17 

4.74 

2.26 

4.75 

1.88 

4.56 

2.08 

OVII(o,574)/(o.6977) 

12.83 

4.96 

3.39 

5.06 

6.66 

6.06 

12.66 

12.26 

12.67 

3.62 

12.39 

2.93 

12.54 

2.90 

12.56 

2.31 

9.41 

5.87 

OVlll(Lya)/(Lyl3) 

4.61 

2.57 

1.37 

1.29 

1.14 

1.38 

5.27 

1.43 

5.14 

1.47 

1.47 

1.82 

5.14 

2.19 

2.40 

2.01 

1.18 

2.57 

OYlll(^Lyo,) / (Lyy) 

12.47 

3.88 

133 

1.50 

1.34 

1.52 

14.35 

2.14 

13.96 

1.64 

2.74 

1.95 

13.86 

2.38 

5.64 

2.16 

1.50 

2.92 

OYlll(Lya)/(LyS) 

25.52 

4.60 

1.36 

1.67 

1.62 

1.64 

29.89 

3.44 

29.02 

1.79 

5.03 

2.07 

29.16 

2.59 

11.38 

2.35 

2.18 

3.18 

MgXI(1.352)/(1.579) 

4.85 

1.89 

2.82 

1.64 

2.42 

1.32 

487 

4.34 

5.36 

1.43 

4.49 

1.72 

4.94 

2.27 

4.88 

2.15 

2.17 

2.59 

MgXI{1.352)/(1.660) 

13.05 

1.95 

6.90 

2.90 

5.66 

2.28 

13.0 

11.36 

14.33 

1.98 

11.84 

2.00 

13.11 

2.61 

12.89 

2.46 

4.75 

3.06 

SiXIII(i,865)7(2.182) 

4.99 

1.62 

4.72 

1.55 

1.63 

1.62 

5.02 

2.82 

5.05 

2.04 

3.26 

2.34 

5.14 

2.99 

4.28 

2.85 

1.64 

2.64 

SiXIII(i,865)7(2.294) 

13.05 

2.28 

12.2 

2.25 

2.88 

1.98 

13.23 

6.64 

13.32 

2.34 

7.96 

2.63 

13.64 

3.47 

11.03 

3.29 

2.75 

3.05 

SiXlV / (Lyy) 

13.10 

7.18 

14.43 

2.06 

2.24 

1.66 

14.41 

2.64 

14.09 

2.11 

5.66 

2.11 

13.84 

2.26 

5.89 

2.18 

3.52 

1.79 

SXV(2,461)7(2.883) 

4.86 

2.20 

5.05 

1.66 

1.76 

1.72 

4.82 

2.36 

4.89 

2.39 

2.96 

2.47 

4.89 

2.90 

3.59 

2.77 

1.83 

2.36 

SXV(2,461)7(3.032) 

13.14 

4.83 

13.7 

2.64 

3.27 

2.21 

13.36 

5.47 

13.51 

2.79 

7.42 

2.84 

13.93 

3.48 

9.66 

3.32 

3.68 

2.81 

FeXVII(o.8211) 7(0.8984) 

- 

4.15 

10.8 

3.9 

5.21 

3.68 

- 

11.5 

- 

4.12 

11.87 

4.12 

- 

4.77 

12.29 

4.32 

5.59 

3.03 

Fexvii(o,8211)7(1.008) 

- 

2.07 

3.20 

1.95 

1.93 

1.90 

- 

3.39 

- 

1.83 

3.43 

1.76 

- 

1.65 

3.51 

1.54 

2.00 

1.18 

FeXXV(6,700)7(7.881) 

4.90 

3.82 

5.20 

4.11 

3.95 

3.4 

5.12 

2.00 

4.71 

2.52 

3.76 

2.4 

3.86 

2.60 

2.77 

2.49 

5.92 

5.33 
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